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AN ELEMENTARY EXPOSITION OF THE THEORY OF THE GAMMA 

FUNCTION. 



By J. L. W. V. 
Authorized translation from the Danish by T. H. Gronwall. 1 

Although there exists a very considerable literature concerning the 
Gamma function, it is but recently that a connected exposition of the 
properties of this function has appeared, with a fairly complete bibliog- 
raphy 2, 3 . However, the monograph referred to does not seem to lay any 
stress on a rigorous exposition, and the treatment is neither carried out 
from any definite point of view, nor based on the simplest foundations. 
The present paper may therefore not be entirely devoid of interest, 
since it gives a short, but rigorous and fairly complete theory of the 
Gamma and allied functions, based on the elementary theory of functions, 
i. e., that part of the theory of infinite series and products, and particu- 
larly power series, which may be treated in a simple and natural way 
without the aid of the calculus. Thus I shall deal below with the proper- 
ties of certain functions of a complex variable, while the connection with 
definite integrals must fall beyond the scope of the present paper; the 
transition to these may however be made with the utmost facility at 
almost any stage, as I may perhaps show on another occasion. 4 

1. Definitions. When we attempt to determine a single- valued function 
/(s) satisfying the functional equation 

(1) /(« + 1) = sf(s), 

it is necessary to inquire first if this problem has more than one solution. 

i Gammafunktionens Theori i elementser Fremstilling. Nyt Tidsskrift for Mathematik, 
Afdeling B, vol. 2 (1891), pp. 33-35, 57-72, 83-84. 

It is hoped that a translation of this monograph, which combines the merits of conciseness 
and elegance in an unusual degree, will be of general interest. The footnotes enclosed in square 
brackets [ ] have been added by the translator, and contain additional references to the literature, 
some propositions supplementary to the text, and finally explanations, in the infrequent cases 
where such have been thought advisable. 

'Brunei, G. Monographie de la fonction Gamma. Memoires de la socigte' des sciences 
physiques et naturelles de Bordeaux, ser. 3, vol. 3 (1886), pp. 1-184. 

* [Of the literature published since the appearance of the original, we shall only mention here 
Brunei's article on definite integrals in Encyclopedie der Math. Wissensch. and particularly 
Nielsen, N., Handbuch der Theorie der Gammafunktion, Leipzig, Teubner, 1906. The latter 
contains a very complete bibliography.] 

4 [This project was unfortunately never carried out.] 
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THEORY OF THE GAMMA FUNCTION. 125 

Let F(s) denote, for the moment, some definite and single-valued solution, 
and write /(s) = p(s) -F(s); it is then seen at once that the relation 

p(« + 1) = p(s) 

constitutes the necessary and sufficient condition that f(s) shall satisfy 
(1). In other words, whenever we have found one solution of the func- 
tional equation, there exists an infinity, since any solution retains its 
character as such upon multiplication by an arbitrary periodic function 
having the additive period unity, and conversely it is evident that every 
possible solution is obtainable in this manner. 

We 6 are led to the discovery of a particular solution of (1) by the fol- 
lowing considerations. 

For any particular solution F(s) that may exist, we have, by (1), 

F(s + 1) F(s + 2) 

F( - s) " i - w+v 

(2) 

F(s + n) 



«(• + 1) • • • (« + n - 1) 



where n, as throughout the following discussion, denotes a positive integer, 
and s a variable which never equals a negative integer or zero but may 
take any other value real or complex. 

Moreover when s is a positive integer, we have 

F(s + n) = (s + n - l)(s + n - 2) • • • n • (n - l)IF(l), 
and therefore 

Therefore, since 

F(s + n) F(s + n) (n - 1) In' 



(&) s(s + 1) • • • (s + n - 1) (n - l)\n' ' s(s + 1) • • • (s + n - 1) ' 
we have, at least for positive integral values of s, 

(c) F(s) = F(l) • lim , , J?~} )] % ^n • 

w w w «_»«, s(s + 1) • • • (s + n - 1) 

But, as we proceed to show, the second factor on the right has a defi- 
nite value for every finite value of s, real or complex, except 0,-1, 
— 2, • • • ; it being understood that here and in what follows n' is the num- 

* The discussion from this point to that similarly marked below replaces the corresponding 
discussion in Jensen's paper. The proof of the uniform convergence of the infinite product which 
defines T(s), however, is that given by Jensen. 
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ber e 8log ", where log n denotes the real value of the logarithm of n. For 

(n -!)!«• ltj^ V 1 +;J 



and the infinite product which is the limit, as w approaches <», of the right 
member of this equation, converges uniformly for every s bounded in 
absolute value and different from 0, — 1, — 2, • • •. An infinite product 

of the form H (1 + a„) will converge if the series ]£ a, is absolutely con- 

r=l ►— 1 

vergent and therefore, in particular, if j» 2 |a„| approaches a finite limit. 
But, given any number N > 1, we have for any s such that \s\ < N and 
any v > N 



('+*)• 



i+; 



- 1 



«(« - 1) 1 s(s - 1)(« 

2! V ' 3! 


-2) 1 


Nl 





and the product of this last expression (which is independent of s) by v 2 
approaches the limit $N(N + 1) as v increases indefinitely. Therefore 
the infinite product under consideration converges, and that uniformly, for 
■every s bounded in absolute value, and this remains true when s is zero 
or a negative integer provided the corresponding infinite factor is omitted 
from the product. 

The function defined by this infinite product is denoted by r(s), so that 



x*( 1+1 J 



<3) roo-jn 

V 

or 

(30 T( S ) - lim -, ,J?~} )l !f tt- 

J _>«, s(s + 1) • • • (s + w - 1) 

It has the properties 6 

6 Weierstrass, K., Ueber die Theorie der analytischen Facultaten. Journ. ftlr Math., vol. 
•51 (1856), pp. 1-60. Reprinted in Abhandlungen aus der Functionenlehre (Berlin, 1886), pp. 
181-262 [and Mathematische Werke, vol. 1 (Berlin, 1894), pp. 153-211]. 
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<4a) T(s + 1) = sT(s), 

from which it also follows that 

r(l) = 1 and Y{n) = (n-1)! 
-For, by (3'), 

T(s + 1) [ (n - Dln^ 1 (n-l)ln' "1 

r(«) SL(s+D(s+2)...(s + n) : s(s + l)"-(s + n-l)J _s ' 

and therefore T($ + 1) = sT(s). But, this having been established, we 
have, as in (2) and (6), 

(d) T(s) = r (* + "> . (n-l)!n> 

w 1W (n-l)ln' s(s + l) ••• (a + n-1)' 

and therefore, since the Umit of the second factor on the right as n ap- 
proaches * is itself r(s), we have, for all values of s, 

r r ( s ± n) 1 
^(n-l)bi'- 1 - 

That r(l) = 1 is seen at once from (3). It also follows from (4a) and 
(4&), since when s is a positive integer we may, by (4a), replace T(s + n) 
by (s + n - l)!r(l) in (46), and (46) then becomes r(l) = 1. 

Conversely the two conditions (4a), (46) serve to define r(s). For (d) 
follows from (4a) and from (d) in turn it follows, by (46), that T(s) is the 
function defined by (3'). 

Thus r(s) is itself a particular solution of the initial functional equation, 
/(s + 1) = sf(s), (1). Notice also that if in (c) we replace F(s) by T(s) and 
therefore F(l) by 1, (c) reduces to (3'). 

Multiplying 5 (3') by 

1 = lim (n + a)' In", 

n— >x 

where a is independent of n but may depend on s, it is seen that 

which is slightly more general in form than (3') and will be used later. 

Before proceeding further, we shall transform (3). It is well known 
that l/(v + 1) < log (1 + 1/v) < 1/v, or 



0< -- log (l +-)<-- -4rT> 

V °\ V ) V V + 1 
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whence it follows that 2(1/V — log (1 + 1/v)) is convergent. On the 
other hand, 

E(;-log(l+;))-l+* + ••• +i-log(n + l), 
and hence we may write 

= lim[l+i+ ••• +i_log(n + l)l, 

which is known as Euler's constant. 7 Its value correct to ten decimals is 
0.5772156649 • • • . 8 Consequently 

r--e* V ^ J J =n(l+- 1 )e^, 
and dividing this by (3), we obtain 

© _^ = e <.. s n(i+5)»-'-. 

As stated above, the condition (4fe) was given by Weierstrass in his famous 
memoir on analytic factorials, 6 which may be considered as an unsurpassed 
model of rigorous and clear exposition. Equation (3) is due to Euler, 9 
who unfortunately soon replaced this excellent definition by definite 
integrals 10 ; in consequence, several of the formal properties of the Gamma 
function escaped his attention. (3') is due to Gauss" who undoubtedly 
was not familiar with Euler's expression. (5) is due to [Schlomilch and] 
Newman 12 , and was also established by Weierstrass. 6 Euler's constant, 
referred to above, is sometimes, without any justification, called Mascher- 
oni's constant. What we have denoted here by T(s) according to 

* Euler, L. De progressionibus harmonicis observationes. Comment. Acad. Petrop. vol. 7 
(1734-1735, published 1740), p. 156. 

» [The arithmetical character of this constant is entirely unknown; it has not even been 
shown to be irrational.] 

9 Euler, L. Letter to Goldbach, Oct. 13, 1729. Correspondence math, et phys. de quelques 
celebres geometres du 18 e siecle, publiee par Fuss, vol. 1 (St. Petersburg, 1843), p. 1. 

10 Euler, L. De progressionibus transcendentibus seu quarum termini generates algebraice 
dari nequeunt. Comment. Acad. Petrop., vol. 5 (1730-31, published 1738), pp. 3(H57. 

11 Gauss, C. F. Disquisitiones generates circa seriem infinitam etc. Comment. Gotting., 
vol. 2 (1813), pp. 1-46. Reprinted in Werke, vol. 3 (1876), pp. 122-162. [German translation 
by H. Simon, Berlin, 1888.] 

° [Schlomilch, O. Einiges fiber die Eulerschen Integrate der zweiten Art. Grunert Archiv, 
vol. 4 (1844), pp. 167-174.] Newman, On r(o) especially when a is negative. Cambridge and 
Dublin math. Journal, vol. 3 (1848), pp. 57-60. 
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Legendre, 13 was written [s — 1] by Euler, U(s — 1) by Gauss and 1/Fc(s) 
by Weierstrass. 

2. Properties. From (5) it follows immediately, by well known ele- 
mentary propositions, that l/r(s) may be expanded in a power series con- 
vergent for all finite values of s, or that l/r(s) is an entire transcendental 
function, with the simple zeros 0, — 1, — 2, • • -. 6 The function T(s) is 
consequently a single-valued analytic function of s with the poles of order 
one 0, — 1, — 2, • • •, and it is seen at the same time that the equation 
r(s) =0 has no roots. All this also follows from the fact that the product 
on the right side of (4) was shown to be uniformly convergent. 

Since sin its has the zeros 0, ± 1, ±2, • • •, one is led to consider the 
function l/(r(s)r(l — s)) which has the same zeros. By reason of (4a) 
and (3) or (5) it is found that 

1 1 °° / s 2 \ 

r(«)r(i - s) = ~sr(s)r(-s) =s ¥v 1 ~^j" 

This product, however, equals (1/x) sin xs, by a well known theorem due 
to Euler and- proved in the elements of the theory of functions, and 
consequently 

(6) r(s)r(i -s) =-^—, 

which is also due to Euler. 14 For s = §, this formula gives 

(7) r(j) = Vx, 

where the square root is positive, since the product expansions (3) and 

(5) both show that when s is real and positive, the same is true of T(s). 
Equation (7) was discovered by Euler as early as 1729. 9 Since it is 
known that 15 

u [Legendre, A. M. Recherches sur diverses sortes d'integrales d^finies, M6m. de Hnstitut, 
vol. 10 (1809), pp. 416-509.] 

14 [Euler, L. Evolutio formulae integralis, etc. Novi Comment. Acad. Petrop. vol. 16 (1771 > 
published 1772), pp. 91-139.] 

15 {For when h is odd, 

. w . 2x . (h - l)rr . . v . ,2x . . h - 1 

sm-rsin-r- •••sin — ,— — = sin J t sin 2 -r- • • • sin 8 „, t 

(a) 

= 2 -/*(l-cos^)(l-cos^)...(l-cos(^). 

But h being odd we have 

(ft-D/2 / 2i — 1 \ 

K* + l = (a+l) n (s 2 -22:008^—^+1] 

and if both members of this identity be divided by z + 1 and x be then set equal to 1 in the result 
we obtain 

(6) h =2«->/* (l -cos£) (l -cos£) ••• (l - cos^^) . 

By substituting (6) in (a) we at once obtain the formula of the text. The proof for h even is 

similar.} 
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oi-i 7. . x . 2t . (h — l)ir 
2 1 * • h = sin -r sin -r- • • • sin - — r , 

where h is a positive integer, we find by making s = 1/h, 2/h, • • • , (h — l)/k 
in (6), multiplying the resulting equations and extracting the square root, 

» r (s) r (s)- r ( ?L T 1 ) =(2 ' )K "°''*-"'' 

which is also due to Euler. In § 4, we shall give another proof of this 

equation without using the properties of the sine. Writing 

( v r(w)w g 

KS,n) ~«(« + l) ■•.(« + n- 1)' 
it follows from (3') that lim (s, n) = T(s); it is readily seen that 

n->oo 

(s, n) ( s + \, n) ...(«.+ *-j^, n) r(n)¥H) ^ 

(/is, /in) ~ T(n/i) " ft ' 

and since the left side of this equation has a definite limit for n -> » , 
the same must be true for the right side, so that 

rwr (. + i)... r (. + »^i) 



T(hs) 



- A; • /r A ', 



where the constant A; may be determined with the aid of (8) by making 
s = 1/h, which gives A; = (27r) (A-1)/2 /i 1/2 . Thus we have proved Gauss's 11 
theorem 16, 17 

(9) r(«)r (s + i V • • r (s + —^ ) - /& 1/2 -*»(2ir)<*- 1 >' 2 r(/is). 

16 [The special case h = 2 is due to Legendre. 1 *] 

17 [This may be generalized as follows: h and A; being positive integers, (9) gives 

r (** + j) * A*<'+w«-ifl(2«-)-(*-»/» n'r(«+j+f). 

whence 

n r ( ks + -f ) = A**.+[*(t-i)i/2-*/2(25r)-[*(*-i)3/« n nris + ^+x). 

Dividing this equation by the one obtained by interchanging h and fc, m and v, it is seen that 

which reduces to (9) for k ■> 1. Equation (9') was obtained from the definite-integral expression 
for log r(s) by Winckler, A., Neue Theoreme zur Lehre von den bestimmten Integralen, Sitzungs- 
ber. Akad. Wien, vol. 21 (1856), pp. 389-426. Compare Nielsen,' pp. 196-198 and 326. The 
remark that (9') is a consequence of (9) appears to be new.] 
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3. Evaluation of some infinite products by means of the Gamma function. 
In the notation introduced above, we have for a positive integer m > n. 



(?)' 



(s, m) _ 

(S ' n)= ('^)('+»TT)-(i^)' 

whence, letting m and n increase indefinitely in a quite arbitrary fashion, 



r(s) - i i im 



m 



Assuming in particular that mjn tends toward a definite limit, we have 
(10) Hm(S) , -D»5(l+i). 

By division of Gamma functions of different variables we obtain from 
either of equations (3) and (5) 

(ID r(«)r(fl = frfi + _i_Vi s —\ 

Writing, according to Mellin, 18 R(s) = ais + a 2 s 2 + • • • + a n s n , and 
denoting the n roots of the equation 

P»(l + i2(l/p)) = p n + a lP »- 1 + ... +a n = 

by pi, P2, —-, Pn, it follows from (5) that 

n*r 

r(a - pis)T(a - p 2 s) • • • T(a - pns) 

=^.( 1+B (i))^(l +B (_i_)) e -«./.. 

We note the special case R (s) — — s n , which gives for a = 1 

1 

r(l - s)T(l - e 2,ri >s)r(l - e 4ni i n s) • • • r(l - e l< - n - V2Mii ln s) 

a formula indicated by Liouville. 19 

18 Mellin, H., Eine Verallgemeinerung der Gleichung r(x)r(l — x) = ir/sin rx. Acta Math., 
vol. 3 (1883), pp. 103-104. Ubei gewisse durch die Gammafunktion ausdnickbare Producte, 
ibid., pp. 322-324. 

" Liouville, J., Note sur la fonctkra gamma de Legendre. Comptes rendus vol. 35 (1852) : 
pp. 317-322, and Journal de math., vol. 17 (1852), pp. 448-453. 
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When B(s) is an even function of s, the convergent product 

$0+*(rh)) 

is expressible by trigonometric and exponential functions alone, on account 
of (6). This result is also easily derived from Euler's product expansion 
of sin tx without introducing the Gamma function. 

4. Expansion of log r(s). Any finite and definite limit lim u n may 

n-^oo 

m 

be rewritten as an infinite convergent series Uh + ]C (w„+i — u v ). This 

remark will now be used to transform equation (J"), by taking logarithms 
on both sides and making a = s. It should be noted that in the following 
any logarithm is defined as that branch of the logarithmic function which 
is real for real and positive values of the variable. When s is represented 
in the complex plane in the familiar way, both sides of the equation 
referred to will then be single- valued as long as we do not make a circuit 
of any of the points 0, — 1, — 2, • • •. Thus, writing 

u n = log T(n) + s log (n + s) — log s — log (s + 1) — • • • 

- log (s + n - 1), 
we obtain 

w„+i - Uy = log v + s log - + log (s + v) 

-[(• + *+*) log * - ^ 1 - l] - [(* + \) log^ - l] 

-[(,+|)log S -±^i-(,-i)log S -±- y ]. 

Now the expression inside the last of the three square brackets is the 
general term of a convergent series, since 

lim (n + §) log S n ^ 1 = Jim (n + \) log ( 1 + ^— J = s, 

and the sum of this series for v = 1, 2, 3, • • •, therefore, in consequence 
of the remark made above, equals s — \ log (1 + s). Furthermore, 

2 {v + |) log 1 converges, since for v > 1, 

(' + » 1 °g( 1 +;)" 1 -I5S-ISS+ •"> 
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and the product of this expression by v 2 is bounded. Consequently 

[s + v + 1 1 
(s + v •+• J) log J- 1 must be convergent, since 2(u y+1 — u r ) 

converges, and we may write 

(13) «(«) -£[(« + * + J) log^qp^-l]. 

It is then seen at once that 

tog r(s) = ui + E (Wv+i - «,) 

= [slog (1 + «) - logs] +[«(«) - (« + i) log 1 ^ + l] 

- «(X) - [• - * log (1 + •)] 

= (s — I) log s — s + 1 — co(l) + «(s). 

It remains, to determine the constant 1 — co(l), which we denote by k' 
for brevity, and to this purpose we observe that lim wis + n) = 0, 

since «(s + n) is 2A« remainder in the convergent series co(s), counted from 
the nth term. Making s = n and s — In respectively, it is found that 

2 log r(n) = (2n - 1) log n — 2n + 2k' + 2a(n), 

log T(2n) = (2rc - J) log 2re - 2n + A' + <o(2w), 

whence by subtraction 

*' + 2«(n) - «(2n) - log r(n) ^""' I/ ' 

_ 2 1 / 2 r(n)n 1 ' 2 

~ 10g K* + D---(* + n-l)' 
and increasing n indefinitely 

Jb' - log [2»»r(J)] = log -VS. 

Hence the final result becomes 

(14) log r(s) - (s - i) log s - s + log ^ + w(s). 

The infinite series «(s) is due to Gudermann. 20 As special cases of (14) 
we note, f or s = 1 and s = f respectively, 

<o(l) = 1 - log -yl2lr and «($) = %{l - log 2). 

» Gudermann, C, Additamentum ad functionis r(a) theoriam. Journal fur Math., vol. 
29 (1845), pp. 209-212. 
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We have seen above that w(s) converges toward zero when s increases 
by positive integral increments toward infinity. Therefore, by (14), 
(s — J) log s — s + log V2ir becomes an approximate expression for 
log r(s) when s increases in the manner indicated. When s is a positive 
integer, we have, in particular, the approximation formula known as 
Stirling's 21 formula. Before proceeding to generalize this result, we shall 
give some applications of (14). We have 

„/ s ,- r(n)n« 

r(«) = ]£ «(« + i)...(, + n-i) 



~ S. s(s + 1) • • • (s + n - 1) ' 
and consequently 

(15) rw - ^Hr. »(» + i)-(« + n-D ■ 

which is due to [Enneper and] Gilbert. 22 Furthermore we may determine 
by means of (14), as shown by Cauchy, 23 the limit k of 

T(ri) h n h - 1 l 2 h hn 

T(nh) = hl ' 2 ( 2 *) h ~ ll2eh " (B) - (nW , 

or A; = /i 1/2 (2ir) A ~ 1 / 2 . Thus Gauss's theorem (9), and (8) as a particular case 
of it, is proved anew. The most elegant proof is however obtained from 
(15). The notation w(s).is due to Catalan; Binet 24 denoted this function 
by fi(s) and Cauchy by ^(s). 

5. Investigation of u(s) and log T(s) for large values of 1 5 1. When | a + 1 1 
> 1, we have 

(a + i) log(l +!)-! = - [( a + 1) - J] log(l -^n) ~ 1 



S2K" + 1) (« + !)'" 



21 Stirling, J. Methodus differentialis sive tractatus de summatione et interpolatione serierum 
infinitarum, London, 1730, p. 135. 

42 [Enneper, A., Tiber die Function n von Gauss mit complexem Argument. Diss. Gottingen 
1856.] Gilbert, Ph., Recherches sur le d6veloppement de la fonction r et but certaines int^grales 
d6finies qui en dependent. M6m. Ac. Belgique, vol. 41 (1873), pp. 1-60. 

23 [Cauchy, A. L., Exercices de math., vol. 2 (Paris, 1827), pp. 91-92, and Nouveaux exerciccs, 
vol. 2 (Paris, 1841), pp. 407-408.] 

24 Binet, J., Memoire sur les integrates Euleriennes et sur leur application a la theorie des 
suites, ainsi qu'a 1'evaluation des fonctions de grands nombres. Journal de l'Ecole Polytechnique, 
cahier 27 (1839), pp. 123-343. 
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Now 2~7~~xrn = ~& * or ^ = 2, 3, 4, • • • , and consequently 



(« + i)log(l+±) 



1 
< 



12S|« + 1|' 12 |« + l|(|a + l|-l)' 
Applying this to the series (13), and observing that 

*|« + F + l|(|« + * + l|-l) 

must be convergent since v 2 /[\ s + v + l\(\$ + v + l\ — 1)] is bounded 
for v increasing, it follows first that w(s) may be expanded in the convergent 
double series (Binet 24 ) 

o»( s )=E, ^ 



(16) 

_ly_i_4.i v L__i_ ... 

" 12 £ (« + m) 2 + 12 A (s + M ) 3 + ' 
and second, that 

1 w(s) I < ^S| s + H(|Uh-D' 

all this under the assumption that |s + l|, |s + 2|, •••, are all greater 
than unity, or that s, when represented as usual in the complex plane, 
lies outside a sequence of circles with the radius 1 and the points — 1, 
— 2,-3, as centers, which we shall assume here for the sake of simplicity. 
(If s is inside or on the circumference of one or two of these circles, it 
is only necessary to omit one or two corresponding terms in the series 
for <o(s).) 

Before proceeding to consider the last inequality for complex values 
of s, we shall make an observation connected with Stirling's formula. 
When s is real and positive, we have 

1" 1 1 °° / 1 1\1 

0< co(s)< i2?( s + ; _ i)( s + „) = 12?V« + " - 1 ~s~+~v ) = l2i ' 

and consequently by (14) 

log r(«) = (s - J) log s - s + log V2^ + ^, < 6 < 1, 

a proposition which is frequently used and is commonly proved by the 
integral calculus. 25 

28 [A more accurate formula is readily obtained by evaluating a lower bound for «(s) in. the 
same manner as the upper bound was found above. From (16) we obtain, s being real and positive, 

/ \ ^ ! 5 1 v l« 1 11 



12 f (s + vf 12 f (s + »)(« + i/ + 1) 12 s + 1 ' 
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In the general case, let (1 + p) be the shortest distance from s = x + yi 
to any of the points — 1, — 2, — 3, • • • ; then p is positive according to 
our hypothesis and | s + v \ > 1 + p, so that 

1 1 1 1 +p 1 



|« + f|(|« + H -1)= \s + v\ 2 t l_ p l« + " 

1 + p 
1 + p 1 



P x 2 + y 2 + 2vx + v 2 ' 



First supposing x positive, and noting that A 2 + B 2 .S $(A + B) 2 when 
A and B are real and positive, the last expression will be less than 

1 + p 1 „. 1 + P _2 1 + P 2 

P |s| 2 + " 2= P (M+") 2 P (|s| +"-D(l« i+") 

and consequently 

(I) l-WK^-'g d.i+.-i -TTTn) 

On the other hand, let — (m + 1) < x S — m, where m is a positive 
integer or zero; then 

SlT+71 5 = § | i + H 2 + £. I s + H 2 + J=l. |s + H 2 ' 

where, f or m = or 1, the first sum on the right side should be omitted. 
In the second sum on the right side we have 

and consequently 

logr(s) - ( 8 -|)logs-«+log-sSJ + ^.|^, 0<9<1. 

Retaining two terms in (16), we may even show that < 8 < 1/2 in the last formula, since 

, w i; l . l • l 



12? ( 8 + v)'^12f (« + !«)» 
and, for * real and > 1, 

whence, making x = s + v, 

1 • / 1 1 \ 1 1 

<" (s) > i2 ? lr+7^i " iTT+i ) - 12 • r+r • 

This proof is believed to be new; for a still closer limitation of 8, obtained by definite integrals, 
see Sonin, M. N., Sur les termes complementaires de la formule sommatoire d'Euler et de eelle de 
Stirling, Comptes rendus, vol. 108 (1889), pp. 725-727, and Annales de l'Ecole Normale, ser. 3, 
vol. 6 (1889), pp. 257-262.] 
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S I s + p\* = (x + m) 2 + f + (x + m + l) 2 + y 2 < y 2 + y 2 ' 

making v = m — y. in the first sum, so that n = m — 1, m — 2, • • ♦, 2, 1, 
and p = m + l+/iin the third sum, so that ft = 1, 2, 3, • • •, we find, 
since rc + m^O, x + m-f-l >0, 



_ = y i <y 

12 ^ f ,. — (t. 4- WW -J_ «2 = ^ 



£i I s + v I 2 ~, (/* - (x + m)Y + y 2 _ ^i m 2 + 2/ : 

^ *-' ..2 _1_ a<2 == *"^ 



<2Z 



Sju 2 + 2/ 2== ,S(m+|2/|) 2 
1 



= S 



£J(lvl+/*-l)(lvl+/0 
1 



„=f +2 | s + H 2 ,tl (m + (x + m + l)) 2 + s/ 2 

< Am 2 + 2/ 2< 1 2/ 1 ' 
so that finally, adding the three sums, 

(id i-wi<T(i7r + ?)' (xs0) - 26 

Since (1 + p)/p is bounded, it is now seen immediately from (I) and (II) 
that w(s) converges uniformly toward zero, when the distance between s and 
the nearest point on the negative real axis increases indefinitely. In par- 
ticular, this condition is satisfied when | s | increases indefinitely while 
the arc of s 27 is constant and different from t. When | s | increases in 
the manner indicated, it therefore follows from (14) that 

( 17 ) Um ,r (S) fir- " J » 

where that branch of the multiple-valued function s*~ 1/2 is taken which is 
real for real and positive values of s. This equation may be used to extend 
the domain of validity of (4b) ; it is readily seen that, when t is any complex 
quantity, 

r(« + t) , 



lim 



r(s)s* 



We have used above the expansion of (a + i) log (1 + l/«) — 1 in 

M [The form and proof of (II) have been slightly changed from the original.] 
n [When s is written in the form s = | s | e ei , 9 is the are of s]. 
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powers of l/(a + 1); the expansion in powers of 1/a would have served 
the same purpose. In fact, for | a \ > 1, 

( „ + i)bi ( 1+ i)_ 1 .gt i ^a^ 

whence we deduce, instead of (16), the expansion 

(18) " W -i 2r(r + l) £(. + *)' 

1A 1 1 f. 1 

~ 12£S(« + M ) 2 12£j(« + M ) 8+ '"' 

vaUd for | s |, | s + 1 |, | s + 2 |, • • • all greater than unity. This formula 
is also due to Binet. 28 Neither (16) nor (18) is adapted to the numerical 
calculation of «(s) (and log T(s) through (14)); expansions suitable for 
this purpose will be derived in § 14. 

6. The function f (s). Definition and properties. It follows from (3) that 

log M = _ log ( l+ .«) + £[ nog(l+ I)_^ 1+ _^)], 

which series is uniformly convergent, since it is derived from uniformly 
convergent infinite products. When s is not equal to zero or a negative 
integer, T(s + t) may be expanded in a Taylor series r(s + t) — T(s) 
+ tT'(s) + • • • for sufficiently small values of \t\. Substituting this 
on the left and expanding both sides in powers of t, we find by comparison 
of the coefficients of t 



r(« 



0„_i + £ [log(l+ i)_^]. 



This function, which will be denoted by i]/(s), may be expressed in the 
following forms, using the series 



and 



*-*l>*K)] 

s VVs + J' + l s + vj' 
(19) 

28 Binet, J., Abstract of paper quoted in u , Comptes rendus, vol. 9 (1839), pp. 39-45. 
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These series are uniformly convergent for all finite values of s; for choos- 
ing an N so large that N > \ s |, it is seen that f or v > N 



v + 1 s + v 



s - 1 



(v + l)(s + v) 



< iV + 1 



(i> + l)(v-N)> 



the last expression being independent of s and being the general term 
of a convergent series. Since 



s+n-1' 
we have 
(190 *»-ttn(lo f .-i- r £ T j^^-j). 

and using = lim [log (n + a) — log n], this may also be written 

n— >eo 

From (19) we furthermore obtain 

(19-) m - m = t (&-, - j^,) = 1 9+ \^ , + r) ■ 



which is somewhat more general in form than (19). 
From either of equations (19) it is seen at once that 

*(D = -c, 

and that 
(20) 



ss + 1 s + n — 1' 

Substituting the last expression in (19'), we find 

(21) lim [*(« + n) - log n] = 0, 

which will be generalized below. Conversely, it is obvious that (19') may 
be derived from (20) and (21), so that these two equations are sufficient to 
define \ff(s). 

The expansions of ^(s) being uniformly convergent, ip(s) is a single- 
valued analytic function of s which has poles of the first order at the points 
s = 0, — 1, — 2, • • •, but is holomorphic for all other finite values of s. 
By (19'") we have, for t = 1 - s, 
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= _i_L._i L_ . _i L_ . 

s^l-s l + s" r 2-s 2 + s" 1 "'"' 

and consequently, by a well-known theorem due to Euler,* 

(22) 4,(8) - *{1 ~ s) - - v cot its, 

which could also be derived from (6) in the same way as (19) was deduced 
from (3). Writing, for brevity, 

r i i 1 1 1 

[s, n] = log n 



8 8 + 1 8 + n — V 

it is readily seen that, if h be any positive integer, 

h[hs, hn] — [s, n] - s + ^, n - • • • - H ^— , n = h log h, 

whence, letting n increase indefinitely, 

(23) H(ha) -*(«) -+( 8 +j^ *( s+ ^r0 - fcl °g*- 

These propositions could all be obtained by logarithmic differentiation 
of the corresponding ones for r(s) (Gauss 6 ). I have preferred the above 
elementary deduction, which is as simple as could be desired. What we 
have denoted here by ^(s), according to Cauchy, was written ¥(s — 1) by 
Gauss," <p(s) and z'(s) by Legendre 13 and X'(s) by Binet. 24 

The series for yf,(s) being uniformly convergent, i/,(s + t) may be 
expanded in ascending powers of t for 1 1 \ < \ s |, | s + 1 |, | s + 2 |, • • •, 
and this may be effected by expanding the individual terms on the right 
side of (19). By Taylor's formula, the comparison of the coefficients of 
t m shows that 

(24) *<-»(«) = (- l)-"w! £ (s+ \ )n+1 , (m S 1). 

These functions are readily seen to be completely defined by the equations 
^ <«>(, + 1) - ^ <->(«) = ( "i|" mI , and Urn if,™(s + n) - 0. 

S n— »« 



•{For 

sin « = *$ n ( 1 — 3 ) • 



Hence 



— r cot xs » — t- log sin ir* = — - + 2 I ri I • r 
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They furthermore satisfy relations analogous to (22) and (23), which we 
shall not, however, stop to prove. 

7. Evaluation of some infinite series and products by means of ^(s). By the 
use of the functions ^(s), lA'(s), • • • which we have introduced, it is easy to 
find the sum of series of very general form. For instance, when F(s) is a 
rational function of the form 

as"- 1 + bs n ~ 2 + • • • + I 
w s" + as"" 1 + j8s"- 2 + • • • + X' 

the sum of the convergent series S[F(s + v) — F(t + v)] is found by 
resolving F(s + v) and F(t + v) into a sum of partial fractions in the 
familiar way, thus decomposing the series into several others which may 
be evaluated individually by means of (19'") and (24). When a = 0, the 
series XF(s •+- v) becomes convergent. In the particular case when F(s) 
is a rational function of s 2 , it is seen incidentally that SF(s + v) may 
be evaluated without the aid of ^(s), ^'(s), • • • by the exclusive use of 
Euler's partial fraction expansions for the trigonometric functions. As 
an example of the application of these remarks we may mention 



(s + v — pi)(s + v — p 2 ) • • • (s + v — p„) 

f(s - pi) 



(Pi — P2)(P1 — P») • • • (Pl — Pn) 

f(s — Ps) 

{Pi — Pl) (P2 — Ps) • • • (P2 — Pn) 
$(S - Pn) 



(p» — Pl)(Pn — P«) ' ' • (Pn — Pn-l) * 

where all of pi, p 2 , • • •, p» are distinct.* 



•{Let 
By a theorem of Euler, 



Hence 

CO 



/(*) - n (a -pi). 

(=1 



2 1//'(p«)=0. 



«.ftx._«u,.fij.._«.i" 2 .2 



v=0 (« + »-«)'•• (S + » -Pn) v r Q <*„ /'(p<)(s +y -pi) 

= 1 i r_J l - 1-^-= s 3 r_j___j— i_j_ 



£-«/(« + •>) 
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It is seen from (5) that 
Furthermore, it follows from (19) that 

00 

e «vv(«) _ e -«-«/«TT c «/i'-p/(»+K)] 
1 

and consequently (Mellin 29 ) 

< 2 « ^'rlrT)-n(i+ s -TT) 

a generalization of (5) and reducing to the latter for s = 1. Making 
t = 1, we obtain 

(26) e+w = « ft ( 1 + j-^ ) e- ,/( * + ". 

8. Investigation of ^(s) for large values of | s |. Writing 

«.- log (•! + •)- J -j^ s + l-! ' wo = logs, 

it follows from (19") that f or a = s 

00 

4>(s) = hm u n = Mo + 23 (w„+i — «,,), 

»->00 

and since 

^ +1 -^ = log(l+ s -^)-^, 
it is seen that, making 

(27) „.( 8 ) =£(_!_ _ log ( 1+ -^-)), 

we have 

(28) *(*) = logs - «*(«). 
When | a + 1 | > 1, then 

i-log(l+i) = (a + 1 1) _ 1 +log(l-^) = |:(l-i)^ I? , 

and consequently 

1 /INI"! 1 

-- lQ g ( 1 +-j|<£| a + 1 |.-= | tt + i| ( | tt + i|_i)> 

*» Mellin, H., Om gammafunktionen. Ofversigt Akad. Stockholm, 1883, no. 5, pp. 3-20. 
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whence we conclude that for | s + 1 |, | s + 2 |, • • • all greater than unity, 

a series which we have already studied in § 5. When | s | increases in 
the manner indicated there, «*(s) -» 0, or 

(29) lim ty(s) - log s] = 0, 

which extends the domain of validity of (21). For |s + l|,|s + 2|, 
all greater than unity we find 

(30) +® -U 1 ~\)£vhr- 

On the other hand, the expansion 

•* - logfl +±) - ±^ZV- A, (I oc | > 1), 
a \ a) 2 v a mi 

gives 

which is valid for | s |, | s + 1 |, • • • all greater than unity. Expansions 
more convenient than (30) and (31) for computing w*(s) numerically, 
will be derived in § 14. 

Equation (29) may be used to determine approximately the roots of 
the transcendental equation ^(s) = 0. According to (19), the imaginary 
part of i(>(x + yi) equals 

which cannot vanish unless y = 0, so that all the roots must be real. 
Furthermore it follows from (19"') that 

00 1 

\p(x) — \f/(x') — (x — x') 2 7 — ; — w / i — c , 
rv ' rv ' v ' o (x + v)(x + v)' 

and every term in the series is positive when x and x' are both in one of 
the intervals (0- • • + «), (- 1 . . -0), (- 2- • • - 1), • • • ; since ^(+ °o) 
= + oo, ^( — » ± 0) = =F oo, it is seen that \}/{x) increases steadily 
from — oo to + * in any one of the intervals in question, so that each 
interval contains one root only. The positive root lies between 1 and 2, 
since ^(1) = — C < 0, ^(2) - 1 - C > 0, and its value 1.46163- • • was 
computed by Legendre. Denoting the negative root in the interval 
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(- n- ■ • -n + 1) by - n + x n , (0 < x n < 1), it follows from (22) that 

7r cot irx„ = ^(n + 1 — x„) = log n + 5„, 

where, according to (29), 5„-»0 for ra-» ». Therefore 

1 , 
x n = - arctan , 



t log n + S„ ' 

whence we conclude that a;„logw-*l. For large values of n, the root 
in the interval (— n- • ■ — n + 1) is thus approximately equal to 

1 



n 



logn' 



a result due to Hermite. 30 Remembering that T(x) never equals zero, 
we see at once that T'(x) = T(x)^(x) = must have the same roots as 
f(x) = 0, which is of importance when inquiring how r(a;) varies for real 
values of x. 

9. Evaluation of f(x) for rational values of x. On account of (20), it is 
sufficient to evaluate 

when pfq is a proper fraction (p and q integers, q > p > 0). Instead of 
this expansion we shall consider the power series 



S(t) = Z (-±-j - — £_ V* 

SVv + 1 P + vqJ ' 



which is convergent for t = 1. By Abel's theorem on power series, 81 

" Hermite, Ch., Sur l'integrale Eulenenne de seconde espece. Journal fur Math., vol. 90 
(1881), pp. 332-338. 

» ["When S a v is convergent, the power series F(t) = 2 a v P approaches the limit 1 a, 
oo 

when t increases through real positive values toward unity." The proof commonly given is this: 

CO oo 

2 a, being convergent, the remainder after n terms r„ -= % a„ -* for n -> », and hence there 

exists a positive quantity A such that | r„ | < A for n = 0, 1, 2, • • • and to any c > but as small 
as we please, there corresponds an integer N such that | r„ | < c f orn > N. Obviously 

and consequently f or < t < 1 

F(t) - 2,0, - Sc(fr - 1) - 5 (T - r„ +1 )(f - 1) 

1 1 

= lim [(r» - r t Xt - 1) + (r, - r,)(i* - 1) + • • • + (r. - r„ +1 )(<» - 1)] 

ft — > 00 

= lim [- n(l - - Mt -t>) r»(<»-» - <»)] + lim f*+,(l - P), 

n-9» n— >oc 
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*(f)--c + itaflw. 



We have, however, f or 1 1 \ < 1, 

Z^fl= - *-« log (1 - ««), 
and the problem thus reduces to the determination, for 1 1 \ < 1, of 



P+1-? 



T^P + vq 



, which consists of q times every gth term in the series 



oe fv 

?,-= -log (1-0. 

1 v 

oo 

But given the power series F(t ) = 2 c»^> it is well known that, for = e 2 **", 



f 
q-l 



Z e-wFite*) = q{c p t* + w+« + <w+ 2 * +•••), 



since 



when 



~, «-» - 1 



1 - 1 
S'-p - i 



= 



p is not divisible by q, but 

£ 

*=0 



when v — pis divisible by q. Hence we find, for 1 1 1 < 1, 
£(0 = - «" log (1 - *«) + Z «""* log (1 - *0 M ) 

»=o 

and since r»+i — * 0, we have 

F(t) - 2a„ = - ir„(i»-» - <«-) = - (1 - t) Jr^ 1 - 1 r„(i«-» - <«-). 

1 1 iV+l 

For <<< 1, it is clear that 1 - t > 0, t" _1 — <" > 0, whence 

1 X+l 

< (1 - t) 2 A • 1+ S ^fl"" 1 - <") 

1 2V+1 

= (1 - <)#4 + €«•»• < NA(X - t) + « 

< 2« for < 1 - t < «/iVA, 
can be made as small as we please by taking t positive and less than, but 



so that 



F(t) - 2 a, 
o 
sufficiently close to, unity. 
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e-* bg \~ + £ r*> log (i - to*) - (t*-< - 1) log (1 - t), 

and letting t increase to unity along the real axis 

* (f ) = ~ ° ~ l0g q + ?! r " P l0g (1 " ^ 

Thus our problem is solved, and it only remains to transform the expres- 
sion obtained. Replacing p by q — p, a new equation results, and adding 
it to the above, we find 

*(J) + *( 1 Y ? ) - - 2C - 2 log* + 2 £ cos^log(l - *»)• 

Since ^(p/g) is real, the imaginary part of the right member vanishes 
identically, so that 

*(f) + f( i f 2 )=-2e-21og S +g C08 2 -22 l„g(2-2cos^), 
which, together with 

gives 

,(f)--C-|cotf-lo gS+ ig C os^,o g (2-2co 8 ^). 

Since 

2T(g — jLt) 2x/i . 2x(o — u)p 2irup 

cos — „ = cos — - and cos — ^ ^ = cos — ^, 

? Q q q ' 

ft ^ 1 

the summation on the right side may be extended from y. = 1 to n = *-?>— 
or g according as q is odd or even. The final result thus becomes 

(32) ^(|)=-C-^cot^-log g + f'cos^log(2-2cos^), 

where q' = q/2 or (q — l)/2 according as q is even or odd, and the accent to 
the right of the summation sign indicates that the term corresponding 
to /u = q/2 in the first case should be divided by 2. The proof given 
here is clearer and much simpler than that of Gauss, 12 to whom is due the 
remarkable theorem expressed by (32). As special cases, we find 
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(!) 



+ c 



2V3 



- 2 log 2 

- 1 log 3 



2 a/3 



7T 

2 



- f log 3 

- 3 log 2 
+ | - 3 log 2 

* /i 2" . . _ , a/5. 3+ a/5 
-2V 1 - ^-|log5 + T log- T - 

, » /1 2" . . - , a/5. 3 + 2V5 
+ 2-V 1 - Vl " * lo g 5 + T log— 2~ 



+ 



_2_ 

a/5 



1 + -5 - f log 5 - -j log 



a/5, 3 + 2V5 



10. The functions P(s) and Q(s). In § 2 we saw that T(s) is a single- 
valued analytic function of s, having poles of the first order at the points 
0, — 1, — 2, • • ■, but holomorphic for all other finite values of s. In 
other words, T(s + t) may be expanded in a power series in t for \t\ 
sufficiently small, unless 3 is zero or a negative integer. In the latter case, 
however, we have for sufficiently small values of 1 1 1 



T(- v + t) = m^-V 1 + a + aj, + 



since 



Urn *r(- v + t) = lim j, w, r(1 , + ,f 

We are thus led to form the infinite series 

(-1)' 



(-1)" 



(t-1) 



P(s)=Z; 



(33) lw_ V^ + »)» 

which is uniformly convergent for all finite s and becomes infinite at the 
same points and in the same manner as r(s). For when N is chosen so 
large that N > | s |, we have, for all v > N, 
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(-1)" 



v\{v + s) 



v\(v-N)' 



which is independent of s and forms the general term of a convergent 
series. Hence it follows that P(s) is a single-valued analytic function 
of s in the entire plane, which can become infinite at the points 0, — 1, 
— 2, • • • only, but is holomorphic for all other finite s. For sufficiently 
small values of 1 1 1, we have 

P(- v + t) - ^p-r 1 + a + cut + • • -. 

It is evident, therefore, that T(s + t) — P(s + t) may be expanded in a 
power series in t for any finite s and f or 1 t 1 sufficiently small, and by a 
well-known theorem, T(s) — P(s) is consequently expansible in a power 
series in s which converges for all finite s and is denoted by Q(s). The 
result of the above investigation which is due to Prym 32 may be expressed 
by writing 

(34) T(s) = P(s) + Q(s), 

where P(s) is defined by (33) while Q(s) is an entire transcendental^ function 
of s. We have 

P{s + l) sns) z. (j , + 1)!( „ + s + 1) i 2- ( „ + 1)!( „ + s + 1) 

1 + r(x + D! e ' 

and P(s) consequently satisfies the functional equation 

(35) P(s + 1) - sP(s) - e~\ 

while Q(s), on account of (34) and (4a), will satisfy the analogous equation 

(36) Q(s + 1) = sQ(s) + e~\ 

When the distance of s from the nearest point of the negative real axis 
increases indefinitely, P(s) approaches zero since, denoting this distance 
by p, we have 

1 ' o "! P P 

Hence, a fortiori, 

, om .. P(n + s) n 

» Prym, F., Zur Theorie der Gammafunktion. Journal fiir Math., vol. 82 (1876), pp. 165- 
172. 

** [That Q(s) is transcendental follows most easily from equation (36) below, for supposing 
Q(s) to be a polynomial of degree m, the polynomial sQ(s), of degree m + 1, would equal the 
polynomial Q(s + 1) — e~ l , of degree m, which is clearly impossible.] 
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and by (34) and (46) 

(38) lim ,^5+£ = 1. 

v ' «-»• (n — 1) In' 

That equations (35) and (37) are sufficient to define P(s) (Prym 32 ), and 
consequently (36) and (38) sufficient to define Q(s), is seen in the following 
way, which will also furnish a new expression for P(s). From (35) we 
find upon division by T(s + 1) 

P(s) P(s + 1) e- 1 



r(«) r(« + l) r(& + 1)' 

whence, replacing s by s + 1, s + 2, • • • , s + n — 1 and adding the result- 
ing equations, 
f39) P(f)_ P(s + n) = A 1 



r(«) r(s + n) YT(s + v) 

or 

P(s + n) = «(« + 1) • • • (s + n - 1) 

( 39 ') r ,, ^ 1 

X 



and since 



\P(s) e 1 E s(s + 1) ... (s + ^_ 1) J, 



.. P( n + g) ,. P(n + g) v ,. (n- l)ln« n ^ , n 
lim =y- — j — r = lim j— - — -^r~, X lim ~, — r~x = 0X1 = 0, 
„_»«, T(n + s) *-». (» - l)!w« „_»„ T(n + s) 

(39) gives for n -* °° 

(40) SS-^t ' 



r(«) * rr(« + ») 

or 

1 



(40') P(s) = e" l S / L n / . n 

v ' w 1 s(s + 1) • • • (s + v — 1) 

which is due to [Legendre and was rediscovered by] Bourguet. 34 Making 
P(s) = T(s) - Q(s) in (39), we find 



T(s + n) r(«) rr(«-M 

or 

Q(s + n) = «(« + 1) • • • (s + n - 1) 

(4io r a 1 "1 
x[QW+^E s(s + 1) .., (s + y _ 1) J. 

M [Legendre, A. M., Exercices de calcul integral, vol. 1 (Paris, 1811), pp. 339-343.] Bourguet, 
L., Developpement en series des integrates Euleriennes. Annates de l'Ecole Normale, ser. 2, 
vol. 10 (1881), pp. 175-233. 
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For s ^ it follows from (35) and (36) that P(l) = 1 - r 1 , Q(l) = e" 1 , 
and for s = 1, P(2) = P(l) - e" 1 = 1 - 2e~\ Q(2) = 2e~S etc. It is 
easy to obtain general expressions for P(n) and Qin) from (39') and (41') ; 
we find f or s -* 

P(n) = (n-l)![l-r'"fU 

( 42 ) . i, " " 

Q(n) = (n-iy^Zh- 

From (36) it is seen that 



-v\' 



Q(0) - lim ™± !Llll = lim Od + «> - Od) . g, (1)f 

which we shall learn to calculate in § 13. Making s = — n in (41') it is 
further seen that 

Q{~n) =^>"Q(0)+e-i 
(43) 

*[_» n(n - 1) ^ n(n - l)(n - 2) + n! J' 

11. Solution of a functional equation. The form of equations (35) and 
(36) furnishes an occasion to consider the more general functional equation 

(44) f(s + 1) = sf(s) - R(s), 

where R(s) is a polynomial in s. Let S(s) be some definite solution, then 
S(s + 1) = sS(s) — R(s), and consequently 

/(* + 1) - S(s + 1) = sLf(s) - S(«)], 

or according to § 1, /(s) — S(s) = p(s)r(s), where p(s) is an arbitrary 
periodic function with the additive period unity. The most general 
function satisfying equation (44) is therefore 

/(«) = S(s) + p(s)T(s). 

To find a particular solution S(s) we divide (44) by T(s + 1), whence 

M _ f(8 + i) _ TO 
r(s) r(8 + l) r(« + l) ' 

and consequently 

/(«) _ /(* + n) _ A P(s + v - 1) 

r(«) r(s + n) r r(« + iO • 

In this equation, we now let n increase indefinitely, the resulting infinite 
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series converging rapidly as shown by (46), and we obtain 

/(«) _ ,. /(« + n) _ ^ Rjs + v-l) 
r(«) ^r( S + n)~V T(s + v) ' 

where the limit on the left must exist and be finite. For the particular 
solution f(s) = S(s) we assume this limit to be zero, 35 or 

(45) is ixm^) if. (^"w = °- 

Hence we find 

W r(«)""r T(s + v) ' 

Any polynomial of degree m may be written in the form 

R(s) — a + ais + a 2 s(a — 1) + • • • + a m s(s — 1) • • • (s — m -f- 1), 

as is readily seen by the method of undetermined coefficients, and since 
we have 

R(8) 00,01,02, , dm 

— — r Tvrr + iv7 r\ + 



r(« + l) " r(« + l) ^ r(«) T r(»-i) T ^ r(« - m + 1) ' 

it follows from (40) that 

^ ISfr + p - 1) A P(s - m) 

£ r(s + .) -«£«* r(« -ju) • 

But in consequence of (39), P(s — /*)/r(s — /z) noay be written 

P(s-n) = PM + e -iY,. 1 



r(« - M ) " r(«) ^ * ^rfc-M**)' 

and it is therefore seen at once that 

fg(8 + y-l) eP(s) 



E«.+E»»£- 



,ti r(s + v) ~ r(a)^"*^;f=i w '\t{r(*-M + *')' 

or by (46), upon multiplication by r(s), 

m m m ii— 1 

(47) S(«) = eP(«) Z a M + Z •. + Z a, £ (« - 1) (« - 2) • • • (s - p + v) , 

jt=0 /i=l (1=2 »=1 

whereby the proposed functional equation is completely solved. Lind- 
hagen 38 has derived this solution in a somewhat more complicated manner 
and without giving it in explicit form. 

* [Carrying out a suggestion by the author, the translation differs here slightly from the 
original,' in which the limit (45) is assumed to be a constant K, not necessarily zero.] 
« Lindhagen, A., Studier ofver Gammafunktionen. Thesis, Upsala, 1887. 
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It would lead us too far here to discuss the solution of the more general 
functional equation 

f(s + 1) = r( S )f(s) + B(«), 

where r(s) and B(s) are rational functions, and regarding which the 
researches of Mellin 37 may be consulted. 

12. Investigation of P(s) for real values of s. It appears immediately on 
inspection of the expansion (40') that P(x) is positive and decreasing for 
real, positive and increasing values of x. When x falls within the interval 
(0- • -1), we therefore have P(x + 1) > P(2) = 1 - 2tr 1 , and by (35) 

n , . eP(x + 1) + 1 e - 1 
eP(x) = x > ~ x —, (0 < x < 1), 

an inequality which will be used below. In general it follows from (39') 
that 

eP{x) + 1 
x - 1 ' 

eP(x) + x 



(a) 


eP(x - 1) = 


(0 


eP(x - 2) = 


(7) 


eP(x - 4) = 



(x-l)(x-2)' 

eP{x) +x + (x-l)(x- 2) 2 
(x-l)(x-2)(x-3)(x-4)" 

For an x in the interval (0- • -1), the numerator on the right in (a) is 
positive, the denominator negative, and hence P(x) is negative in the 
interval (— 1- • -0); on the other hand, both numerator and denominator 
are positive in (0), and consequently P(x) is positive in the interval ( — 2 • • • 
— 1). To show that the same is true for (y), we use the inequality given 
above, and obtain, the denominator being positive, 

eP(x -A)> e-l + x> + x(x-l)(x-2y 
ef[x ±) > x(x _ 1)(x _ 2)(a . _ 3) (a . _ 4) . 

Making x = 1 — y, < y < 1, the numerator of the last fraction equals 
e - 1 + (1 - yY - 2/(1 - y)(l + y)> s e - 1 + (1 - yY - \{1 + vY 

= e - \y - f 2/(1 - y) 
^e-{-&=e-\\>% 

"Mellin, H., Zur Theorie der Gammafunktion. Acta Math., vol. 8 (1886), pp. 37-80. 
ljber einen Zusammenhang zwischen gewissen linearen Differential- und Differenzengleichungen, 
ibid., vol. 9 (1887), pp. 137-166. [Zur Theorie der linearen Difierenzengleichungen erster 
Ordnung, ibid., vol. 15 (1891), pp. 317-384. tjber den Zusammenhang zwischen den linearen 
Differential- und Differenzengleichungen, ibid., vol. 25 (1902), pp. 139-164.] 
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since y(l — y) = | — (§ — y) 2 s= J. Hence P(x) is likewise positive in 
the interval (— 4- • • — 3). From (/?) it is seen immediately that when 
P(x) is negative for a negative value of x, so is also P(x — 2), and conse- 
quently P(x) is negative in the intervals (— 3- • • — 2), (— 5- • • — 4), 
...(_ 2n - I--- - 2n)---. 

It only remains to investigate the intervals (— 6- • • — 5), (— 8- • • 

— 7), • ■ ■. When <= > but as small as we please, (33) gives 

P (- 2n + *=W)U + ---'> f(-2n + l- e )= (2n _ 1 1) , e +---. 

P(x) must therefore be positive in the latter intervals, provided x is sufficiently 
near either end of the interval. But we may show that P(x) will also assume 
negative values in these intervals. We have by (7) and (/3) 

en +>' f-M-V- "i-f •§•*•*■¥' 

and since, by (40 ), 

«P(J) =2(l+j + F ^ I +..-)<2e, 

P( — V") i s visibly negative, and hence also P(— J ^-), • ■ ■ P( — 2n + §), 
• • • . That the negative values which P(x) thus assumes in the intervals 
(_ 6- • • — 5), (— 8- • • — 7), • • •, cannot be large, appears from the 
equation 

v ' ex x ' 

which shows that when P(x), P(x + 1) and x are negative, we must have 
P(x) - 1/ex > or I P(x) \ < 1/e | x | . 

Since P(x) is a continuous function of x in all the intervals in question, 
it is clear that the equation P(x) =0 must have at least one root in each 
of the intervals (-*£•■• -5), (-6----^), ••• (-2n + §--- 

— 2n + 1), (— 2ra- • • — 2n + i), - • • • but none in the remaining intervals 
(Bourguet 38 ). 

Note. It is a defect of the preceding investigation that it gives no 
information regarding the possibility of more than one root in the interval 
considered. 39 Bourguet 40 has attempted to show that the equation 

88 Bourguet', L., Sur la fonction Eulerienne, Comptes rendus, vol. 96 (1883), pp. 1307-1310, 
and Acta Math., vol. 2 (1883), pp. 296-298. 

39 [That each of the intervals in question contains exactly one root has been shown recently 
by C. N. Haskins, On the zeros of the function P(x), complementary to the incomplete Gamma 
function, Trans. Am. Math. Soc, vol. 16 (1915), pp. 405-412. For a much simpler proof, which 
also shows that there are exactly four imaginary roots, see T. H. Gronwall, Sur les zeros des 
fonctions P(z) et Q(z) associees a la fonction gamma, Ann. de 1'Ecole Normale, 1916.] 

40 Bourguet, L., Sur la theorie des integrates Euleriennes, Comptes rendus, vol. 96 (1883), 
pp. 1487-1490. 
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P(x) = cannot have more than 4 imaginary roots, but his proof is not 
correct. We may note in this connection that it is still undecided whether 
the equation Q(s) = has any roots or not. 41,42 That there are none 
when the real part of s — 1 is negative, was shown by Lindhagen. 36 We 
cannot enter further upon these questions which require more powerful 
methods than those at our disposal here, but on another occasion, the 
author proposes to examine the equation P(s) = k for k real and show 
that for k = 0, the equation has exactly 4 imaginary roots. 

13. Calculation of the coefficients in some power series. By (19) we have 

w + . ) __c + £( r L T -_l- rt ), 

whence, expanding each term in powers of s for \ s \ < 1, 

(48) f (1 + s) = - Si + S 2 s - &s 2 + S t s 3 

(480 Si = C, S n = l+~- n + ^ + •••, (n>l). 

The coefficients 5 2 to jS 35 were computed by Legendre 43 to 16 decimals 
(15 correct) and 5 2 to S 70 by Stieltjes 44 to 32 decimals (30 correct). To 
determine the coefficients in the expansion 

(49) r(l + s) = 1 + ais + a 2 s 2 + ••-, 

which will be convergent for \ s \ < 1, we recollect that 

I"(l + s) = *(1 + s)r(l + s), 
or 

oi + 2a 2 s + 3a 3 s 2 -\ = (- Si + S 2 s - &s 2 H )(1 + a t s + a 2 s 2 H ) 

whence, comparing the coefficients of s" on both sides, 

(49Q (n + l)a„+i = - S t a n + S*a n -i -•••+(- l^+^+i, 

41 Later note. In a paper read before the Mathematical Society of Copenhagen in 1893, the 
author has given approximate expressions for the roots of Q(s) «■ 0. This investigation has not 
yet been published [but will appear in an early number of these Annals]. 

42 [The existence of an infinity of zeros of Q(s) may be proved in the following manner, which 
is not elementary, but has the advantage of brevity. Assuming Q(s) =0 to have only a finite 
number of roots, «i, ss, • • •, s» (or none), the only roots of Q(s) = e -1 will be 1, «i + 1, s» + 1 
•••,*» + 1 since, by (36') , Q(s) = e -1 is equivalent to (s — l)Q(s — 1) =0. But by a theorem 
due to Picard, an entire function which takes each of two different values (here and e _1 ) only a 
finite number of times, is a polynomial. Since Q(s) cannot be a polynomial,' 3 it follows that the 
equation Q(s) =0 has an infinity of roots.] 

43 Legendre, A. M., Traits des fonctions elliptiques et des int^grales Eul6riennes, vol. 2 (Paris, 
1826), p. 432. 

44 Stieltjes, T. J., Table des valeurs des sommes S k , Acta Math., vol. 10 (1887), pp. 299-302. 
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which may be used to compute the coefficients a„ numerically. For 
n = 1, 2, • ■ • we find successively 

a x - - Si, a 2 = KSi 2 + St), <h = - i(Si 3 + ZS& + 2S 3 ), ■•-, 

which however are less suited to numerical computation than (49'). 
Writing (— l)"a„' instead of a n in (49'), the alternating signs disappear, 
and it is seen that if a/, a 2 ', • • • , a„' are positive, this must also be true 
of a'n+i. But ai' = Si > 0, hence all a n ' are positive, and a n has the 
sign of (— 1)*. It is also easy to calculate the coefficients of 

(50) P(l + «) = b + b lS + b 2 s* + •••, (|«| < 1), 
by expanding each term of 

p q+ , >-W+i+.)- 

We obtain 

(500 i.-1-r-, i-f ^ , 

and these expressions are quite convenient for numerical computation. 
The infinite series for 6„ consists of numerically decreasing terms with 
alternating signs, and consequently b n has the same sign as the first term 
of the series, or the sign of (— 1)". 

Since Q(l + s) = r(l + s) — P(l + s), the coefficients in the power 
series 

(51) Q(l+«) = Co + c lS + c 2 s 2 + ..., 

which converges for all finite s, are given by 
(51') Co = e -1 , c n = a n — &„. 

Since Ci = Q'(l), this also determines Q(0), as we saw at the end of § 10. 
The calculation of the coefficients in the expansion 

(52) f(TTl) = ! + *« + *>* + •••» 

which will converge for all finite s, may be effected by comparison of the 
coefficients of s n in 

(1 + a t s + a z s 2 + ■ ■ -)(1 + dis + d** 2 +•■■) = 1, 
whence 
(52') d n + d„_iai + d*_ 2 a 2 + • • • + o„ •» 0, 

which gives, using the previous expressions for a%, 02,03, • • - , 

d t = St, d 2 = %(Si 2 - Si), d 3 = US1 3 - aSA + 28,), • • •. 
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Having once calculated the coefficients oi, a 2 , a 3 , • • • , (52') serves to deter- 
mine di, d 2 , ck, • • •• When an independent determination is preferred, 
we may proceed as follows. By Taylor's theorem and (52) we find 



r(i + « + t) r(i + s) + *r'(i +«) + ... 

= l + d 1 (s + 1) + d 2 (s + ty + ■ ■ • 

and comparing the coefficients of t on both sides of the last equality sign, 

i r'(i + s) l 

" roT7) • r(i + S ) = _ r(i+8)- * (1 + s) = * + 24« + S4a , + • • •, 

or 

(1 + dis + d 2 s 2 H )(Si - S 2 s + S 3 s 2 -•••)= di + 2d 2 s + 3d 3 s 2 H , 

whence finally 

(52") (n + l)d„ +1 = Sid. - S 2 d B _! +...+(- l)»S n+1 , 

which is of the same form as (49') except that the signs of Si, S 2 , • • • are 
changed. The coefficients of the power series for r(l + s) were first 
computed by Jeffery. 46 The coefficients in the expansions of l/r(2 + s) 

= 1/(1 + 8)r(i + s), i/r(«), e<*/r(2 + s), r(2 + s) - (l + «)r(i + «), 

Q(s) and eP(s)/r(s) were computed by Bourguet 34 - 46 to 16 decimals, of 
which the last is doubtful. 

In the papers quoted, the general recurrent formulas are derived in a 
somewhat more complicated manner than above. 

14. Expansions in series of factorials. In the following, we shall expand 
various of the functions introduced previously in series of the form 



s(s + l) ...(« + „- 1)' 
To simplify, we introduce the notation 

(s)„ = s(s + 1) • • • (s + n - 1) 

when n is a positive integer, while f or n = we write (s) = 1. Further- 
more we agree not to let s equal a negative integer or zero. It is then 
seen that 

(S% (S') P+1 _(S%( S' + V\ (aQ 

(s)„ (s) y+1 -'(s) v V s + V ) ~ (s) v+1 {s s) > 



46 Jeffery, H. M., On the derivatives of the Gamma function, Quarterly Journal of Math., 
vol. 6 (1864), pp. 82-108. 

« Bourguet, L., Sur les integrates Euldriennes et quelques autrea fonctions uniformes, Acta 
Math., vol. 2 (1883), pp. 261-295. 
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whence for | s — s' \ > 0, making v = 0, 1, • • •, w — 1 and adding, 

Since (s)» = T(n + s)/T(s), and consequently 

(s')„ r(s) T(s' + n) (n-l)\n> 



n' 



(s) n r(s') (n-l)!w*' T(s + n) 
it is seen at once that, denoting the real part of s by $R(s), 



Um (f0» 





indeterminate, as 91 (s — s') 

00 



> 
= 0. 
< 



On account of (a) we therefore have the expansion, convergent for 
9t(« - s') > 0, 

w «-«' Y(s), +1 ~s i "s(s + l) i 's(s + l)(s + 2)" i " "• 

which is due to Stirling. 21 For 9l(s — s') = the series will oscillate 
(between finite limits) and diverge for 3£(s — s') < 0. We have excluded 
above the case s' = s, where it is seen at once that the series reduces to 
2 l/(s + v) and consequently diverges. 

It does not appear from the preceding that the series is uniformly 
convergent, but this will be shown presently to be true for all values 
of s and s' subject to the conditions: s' finite, 9t(s — s') s € , | s |, | s + 1 |, 
| s + 2 |, • • • s p, where e and p are positive, but as small as we please. 

We suppose these conditions fulfilled and write s = x + yi, 
s' = x' + y'i. The positive quantity A may then be chosen so large 
that A > | $' |, and consequently x' > —A, \y' \ < A. Now taking an 
integer N satisfying the inequality 

N>-A 2 + A, 

€ 

we have f or v s j\T 

v + x'^N+x' >A+x'>0 
and 

\v + S ' | - (v+X f ) <|, 

since 



fr+»o(V i +(^) , -0 < 2-^< 



€ 



- e A* + 2(A + x') 2 
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Therefore 

|y + s|-|" + s'l>!> 

since \v + s\^v + x^v + x' + e. We now write 



u n = 



(«) 



n+l 



Mr, 



M 



n+l 



s + n + 1 



s' + n 



and the preceding inequality becomes, making v = n + 1, 



I «' + n | 



Ur, 



Un+1 



or 



m 



| u n+1 1 < - (| s' + n 

whence 

2 

| W»+l | + | «n+2 | + • • • < - | «' + n 



- I «' + n + 1 | > g , 

- \s' + n + l\-\u n+1 \), 



u n 



9 " 

« 



s' + y 



+ * 



But we have for all values of v, and in particular for v < N, 



and for v S 2V, 



s' + J- 



s + J» 





«' + V 


K A +v 
p 




s + v 


< 


\s' + »\ . 



Consequently (if n > N), 






»' + v 



+ v 



<!n 



?r-f09( 



s' + H + 2 
y + 4 



y + A 
r + A+2 



v + A + 2 



< 



2A(A + 1) -"{A +N - 1) 






i + 



2(* + A) 



) ' 



which is independent of s and s' and converges to zero for n increasing. 
It is thus shown that S«„ or the right side of (/3) is absolutely and 
uniformly convergent when s' is finite, 9J(s — s') a e > and \s\, 

| « + 1 1, | « + 2 1, • • • a P > o. 

On the basis of this result it would be easy to show that the boundaries 
of the regions of conditional as well as absolute convergence of a series of 
the form S cj(s) v are straight lines perpendicular to the real axis, and it 
is also possible to state the necessary and sufficient conditions for the 
expansibility of a function in such a series. As we shall have no occasion 
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to use these results in the present paper, I shall not go beyond this indica- 
tion. 47 

Assuming | « | > \ a' \ and 31 (s — s') > 0, (0) may also be written 

03') ±£i = ± (S% 



Ts' +1 Y(s)„+i' 

The right hand series being uniformly convergent, it is permissible to 
expand in powers of s' and compare with the left side. Writing 

(«')» - CSs'» + Cls'"- 1 + • • • + Ct-xs', CI = 1, 

where the first coefficients C" are readily calculated, 48 we find 

(t) 7h5 " 21 



which will converge absolutely for 9?(s) > 0.* 

47 [An excellent account of this general theory is given by E. Landau, Uber die Grundlagen 
der Theorie der Fakultatenreihen, Sitzungsber. Akad. MUnchen, vol. 36 (1906), pp. 151-218.] 

41 [Using, for instance, the identity («')„ = (a' + n — l)(s')n-i, which gives the recurrent 
formula 

C - ci' 1 + (» - 1)^;!-] 

* {We will give the following direct proof that (7) converges fdr 3J(s) > 0. The ratio of con- 
secutive terms in the series is 

C* m + (» - 1) + c^-vc'tr 1 



C^ _i (M + S + ») M + (« + !-) 

But, as we shall prove, lim C' +M-, /C'"t5 -1 = 0. Hence the series will converge (and that abso- 

lutely) if JR(s) + v - (v - 1) > 1, that is, if 9t(s) > 0. 

That lim C+^/CpH' 1 - may be shown as follows: C^* -1 and C't/f -1 are the sums 

of all the products of the numbers 1, 2, • • •, v + it — 2, taken m at a time, and m — 1 at a time, 
respectively. Hence C£ + ' 4-, may be written 

(w+m-1 „ ZniJh ... nih l <m < «» < • • • < n„< »> +ju - 2. 
To each term «i«» • • • n M in this sum there corresponds a set of n products in C't.'i -1 whose sum is 



«i«j •••»/.( —H- — H 1 - —). 



But if we form this sum for each term in Sni«j • • • n» in turn, each product in C£t£ -1 will be 
taken v— 1 times, since between the number of products in C* +M-1 and that in C^i. 1 } -1 the 
relation pC + JJ -2 ) = (» — 1)(""£!!7 4 ) exists, where (") denotes the number of combinations of n 
letters taken r at a time. We therefore have 

H SWiWj • ' • Tlfi 1 



*\«i wj »„/ » — 1 v v + ii — 2 

But the limit, as p approaches » , of the last member of this inequality is 0. The same must 
therefore be true of C^-VC^tS -1 .} 
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When we have a series 



a v 



a 

converging absolutely for | s | > r, this may be transformed in a series of 

factorials with the aid of (7). We then obtain a double series, absolutely 

convergent for 9fc(s) > r, which becomes, by combining the coefficients 

(Stirling 21 ), 

,. N ,, x a . Oi . % + 01 Oj + 3a 2 + 2ai 

w ^) = (i); + (ir 2 + -(i)r + — w* — + " ,> 

where the general term is 

_A_ = a " + cfa-* + c ^ a -' + — b C:Liai 
( s )i-+i ( s )i-+i 

or symbolically 

0(0 + 1) ♦ • ♦ (o + y - 1) 

(s), + i 

which could also be derived directly from (/3') by symbolic calculation. 
A direct proof that (8) is absolutely and uniformly convergent for 

9t(«) S r + e, 

where € is positive but arbitrarily small, is reached by observing that, 
since 2 a r ls v+1 converges for j s | = r + e/2, there exists a finite and positive 
M such that | a r \ < M(r + e/2)" for v = 0, 1, 2, • • •. Hence 



00 



»+l 



<M V 2A 



00 



H-l 



and the theorem follows by comparison with the series (/3) in which we 
make s' = r + e/2. 
The identity 

1 1 m 

(s + v) m (s + v + l) m ~ (s + v) m+ i 
gives 

' £j (S + »0m+l W L (s)« (S + »)m J ' 

whence we find f or m ^ 1 and n -> 00 (Stirhng 21 ) 

^ 1 11 



00 



„=o (s + f) ro+ i m(s) m * 



It is readily seen that this series is absolutely and uniformly convergent 
for 9i(s) S — A and \s\, | s + 1 |, | s + 2 |, • • • ^ p, where J. is finite 
and p > but as small as we please. 
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Applying the transformation (f) to the series (5), it is seen that, 
making a = 0, 

m +/(s + 1) + ... +„,+„_ 1, . f^.^j. 

"We have, however, 

A v ( s\ s (s + l)_i A v 



v(s + ' 



s + n(s + n + l)„_i (s) r+ i ' 



here 9*(s) > r, s is finite, | s + y. |/| s + n + M | < 1 for n = 1, 2, 
that 

(s + Dv-i 
(s + n + 1)^ ^ ij 
and consequently 



so 



"(s + n)„ 



< 



1 + 



\s + n\ 



A v 



(«) 



H-l 



Since 2 | A„/(s) r +i | is convergent, | 1 + sfv | finite, and | s |/| s + n | -» 
asd-> oo , we therefore find 

M J 

lim £ , " = 0, 

whence finally 

<*) /(«) +/(« + 1) +/(« + 2) + £ -^- 

for 9t(s) > r, which is also due to Stirling. 21 - 49 

These formulas may be applied immediately to the calculation of 
&>(s). It will be remembered from §§4 and 5 that 

«(•) = t[ (« + v + J) log(l +^) - l], 



and that for I s I > 1 



(8 + i)Iog ( 1+ I)_ 1= |; 2l _ttt 

consequently, by (5) and (e), 

(•+*)log(l+J) 



+x 

+ 2)s" +1 ' 



00 t 



roo 



H-l 



rtJ — 



(-D 



H-U 



2(p + !)(* + 2) ^ 2Kv + l) 



+ ( -y:- 1) c ; + 



+ -hCUu 



a [The series to the right in (t?) is absolutely and uniformly convergent for $R(s) |g r + e. 
Tor as shown above, \ A v \ < M(r + «/2)„, whence 



A v 



"(*)» 



<M 



+j+'-i('+SL 



(i^<^(^2 i + 1 )S 



('+iL 



s)„ 



and our statement is proved by comparison with the series 08) in which we make s' = r + «/2. 
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and from (#) we finally obtain 



(53) 



«(«) = Z 



fC v 



7* v(s) v 



which converges absolutely for 9t(s) > 1 B0 and is due to Binet. 24 It is 
evident that \k v \ < T V(1)» = ?!/12. 

60 [The above analysis does not give the complete domain of convergence. Actually (53) 
converges absolutely and uniformly for SR(s) > « > (and consequently equals the holomorphic 
function «(s) not only for 9J(s) > 1 but for SR(«) > 0), but is divergent when 8J(s) < 0. For 
observing that /*/20u + 1)0* + 2) = 1/0* + 2) - 1/2(m + 1), we have 

*■ = -j 1 ( 1 ^ 2+ fe+^) c:: - - -r 1 u°^+is»w->ds, 

or by reference to the identity denning the coefficients C/_ M , 

~ kv= *f<> ( s + i) (*)»*• 
Replacing s by — « in the integral, we obtain, for j> > 2, 

-ft* -J] 1 (J -•)«(! -s)(2 -•)...(!. - 1 - s)efc. 

ButforO < s < 1, we have | 1/2 — * | < 1, « < 1, 1 — * < 1, 2 - s < 2, •••,» — 1— «<»>-l' 
so that 



I J^l-2- •■■(r-l)<to-fr-l)! I 
(■> - 1)(1)k-» _ 1 (l)„-2 



e 



and for SR(s) =2 e and v > 2, 

Kij; "■ h s 1 1 (« + i)„_i r ^ t ' i (s + i)„_i i ' 

in consequence of what was proved above for the series 08), it follows that (53) is absolutely and 
uniformly convergent for SR(s + 1) ^ 1 + « or 5R(s) S «• Now transform the expression for 
— k v by decomposing the integral in two with the limits 0, i and \, 1 respectively, and replacing 

s by 1 — s in the latter. Thus 

J»l/2 
o (J - s)s(l - «)[(2 -«)...(„- l _«)_ (i +,)...(,_ 2 + «)]<&, 

where the integrand is evidently positive, and therefore 

I k * I > X* (i ~ s)s(1 ~ s)[(2 -«)••*("- 1 ~ s) ~ (1 + s) •••("- 2 +«)]<*» 

for a positive 5 less than $ and as small as we please. But f or < s < S < \ we have 1 + * 
<f(2-«), 2+s<3-«, •••, v-2+s<v-l-s, so that (1 +s) •••(»- 2 + s) 
< |(2 — s) • • • (i» — 1 — s), and consequently 



Ky 



> I f! » - *W X - *) ( 2 " •) ' ' ' (" ~ 2 - •) > "~*""" 8 <fe 



Hence 

81(8) < 



tCy 



-j&»(X -»)&-») 



(„ _ 2 -«)-£|(- «)_.,! 



120 



(~ «)>-! 



(8), 



•oo as » -» » when 5R(s) < — «, so that (53) diverges when 
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The first coefficients are ki = J%, & 2 = 0, k 3 = — x%n> *« — ~ *V> 
When s is large, the computation of only a few terms in (53) gives a very- 
close approximation. 

In the same way we obtain from 

i + _j _ i m 2: 2(v + 1)(v + 2)(s + 1)v+1 , 

which is vahd for | s + 1 | > 1, 

hv = 2(p + 1)(f + 2) + 27(7+1) C ' + ' • • + tVC'-i> 
and 

(54) «(«) = £ /S " 



!'»(• + IV 

which is also due to Binet 24 and converges absolutely and uniformly when 

1R(s) ge> 0." The first coefficients are fci = ^, k t — £, A* = -^<nr, &< = ff 

• •. Making s = n, a positive integer, in (53) and (54), we obtain [se; 

(14)] the convergent series 

log (n!) — (ra + 4) log n + n — log V27 

1 1 1 

~ 12n 360n(w + l)(n + 2) 120»(w + 1) (n + 2) (n + 3) 

~ 12(n + 1) + 12(n + l)(n + 2) ^ 360(n + l)(n + 2)(n + 3) 

. ?9 

~*~ 60(n + l)(n + 2)(n + 3)(n + 4) "*" 
Remembering that (§ 8) 

"* (8) -?[n^- l08 ( 1+ s -Tr)]' 

and that for | s | > 1, on account of (5) and («), 

1 i,urfi i.M r ( ~ 1)y+1 -T *- 

- - log ^1 + F j - £ (r + iy+i - £^, 

where 

(_ i)"+i (—1)" 

we find by (#) 

(55) «*(«) = £ 



Ky 



rK«v 



u [Using the method of M , the series is readily seen to diverge for SR(s) < 0.] 
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which will converge absolutely for 9J(s) > l, 62 and in which the first 
coefficients are fa = §, fa = \, fa = \, fa = £# , • • • . Similarly, using the 
expansion 

fVy 



J-^(i+j)-t ( , + 1)( ; + 1) H*-£; 



s W *V ^«/"i (" + !)(« + !) 1 "*- 1 r(s + l)v+i' 
which converges for 9t(s) > 0, we find 

v . v — 1 



fa — _. i i *r , Ci + • • • t" ^CJ_i 



and 

(56) «*(«) = i: * 



which is absolutely and uniformly convergent for $R(s) S e > 0. 68 The 
first coefficients are fa = \, fa = %, fa = *£■, fa = *£$■, • • •. As special 
cases of (55) and (56), we note the convergent series 

l+i + *+ •■• +^zri~logn - C 
1 1 1 



2n 12n(n + 1) 12n(» + l)(n + 2) 

19 



120n(» + l)(n + 2)(n + 3) 
5 



2(n + 1) 12(n + l)(n + 2) 4(n + l)(n + 2)(n + 3) 

469 

120(n + l)(n + 2)(n + 3)(n + 4) '"* 

We shall finally give an application to the function \(/(s). By (19"') we 
have 

W - W8 - s0 =t(_J r: ^-_L-), 

and since for 3J(s — s') > 

_J l = j-. (gO, 

8 - *' S i (s)„+i ' 



K [Actually the series converges absolutely and uniformly for 9i(s) 2: <• > 0, but diverges 
for SR(s) < 0. For we have, by the method of M , 

k v = f* s(l - s)(2 - s) • • • (x - 1 - s)<fe, 

whence < k, < (y — 1) !, and f or < S < 1, 

7 >/.Vl ~ «)(2 - S) ••• fr - 2 - S) "~l~ S ds >||(- *),-i | .] 

» [But divergent for SR(«) < 0, since A„ > J + JCf + • • • + \C v r ^ = §(" - 1)1] 
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it follows from (#) that 
or replacing s' by — s', 

tKT\ ,/ , a ,,s «' S'(S'~ 1) , S'(S' -1)(S' -2) 

which converges absolutely and uniformly for 9J(s + $') S e > 0. 54 As- 
suming s' to be a positive integer, we obtain the finite series 

1 , _1 , , 1 _ n _ n(n - 1) 

s" l "s + l" t "'"" t "s-|-n-l s 2s(s + 1) 

IS. Numerical computation of the gamma function, 55 Let us now consider 
a practical method of computing r(s) for real values of s. It is easily 
shown that r(s) is known when its values in any interval of length J are 
known. 

From (5) we have 



(58) 



L(s) -log TOO = -Cs-log S + E{^-log(l + ^)j. 



It is seen that it is permissible to differentiate the series on the right 
termwise; so that 

w s i [ n s + n J 

This may be written 

L'( S ) = _c+x;j---— ^ -J. 

w Tin s + n — 1 J 

These relations hold for s 4= 0, — 1, — 2, ■ • ■. In general we have 

(59) L"(s) = (-iy(r-l)l± {s + l_ iy , r>l. 

In particular 

(60) 2/(1) = - C, 

(61) LM(l) = ( _ i)'(r-i)! 2 I = (- i)r (r - i)LS r , 

/ft 

where S r has the meaning of (48'). 

Using Taylor's development about s = 1, we have 

log r(s) = L(s) = L(i) + ^--L'd) + ^^-Va) + • • • ; 



64 [But evidently diverges for 9l(s + s') < 0.] 
» This section has been added in the translation. 
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or setting s = 1 + x and taking account of (60) and (61) 

(62) log r(l + *) = - Cx + Z ^- S n x\ 

We show now that this relation is valid for — f == x == J, by proving 
that 

ie m = ^|L<'»>(i + to), o<0<i 

converges to 0, as m -» <» . 
For, if - | =£ a; == §, then 

■ p I < II g "if* xro 1 Iff, 1 11 
|/C "' ~ 1 1 1 + 6x + Y |w + 6x\ m J ' w < 1 1 + Y2-(n - l) m J m~* °* 

As a matter of fact the region of convergence could be shown to be 
— 1 < x 2= 1, but for our purposes it suffices to know that the series con- 
verges f or — J =s a; == §. We have for any a: in ( — §, §) 

log (1 + *)-*-]£(- 1)"£. 



log r(l + x) - - log (1 + as) + (1 - o* + £ (- 1)«(5. - 1)£. 



This on adding and subtracting from (62) gives 

log r(l + x) - - log (1 + 
Changing here a; into — x gives 

log r(l - x) = - log (1 - x) - (1 - Ox + Z 0S„ - 1) ~. 
Subtracting this from the foregoing gives 
log T(l + x) - log T(l - x) 

- " log r^l + 2(1 - C) * - ?2^Ti ^ " 1} - 

From (6), substituting therein r(l + s)/s for T(s) and setting s = x, 



log T(l -J- ar) + log r(l - as) - log ; 



7TX 



Sill 7TX 

This with the preceding relation gives 
log T(l + x) 

= (1 - C)x + i log -j— - - f log . -- 5 I(flWi-l)j 



sin tx 2 6 1 - x 2 V v "^ ' 2m + 1 ' 

This series is due to Legendre. It converges rapidly for == x 2= \ and 
enables one to compute T(s) in the interval 1 == s == f . The value of r(s) 
for other real values of s may then be obtained as already observed. 



